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ABSTRACT 

The observed abundance of giant arcs produced by galaxy cluster lenses and the mea- 
sured Einstein radii have presented a source of tension for ACDM, particularly at low 
redshifts (z ~ 0.2). Previous cosmological tests for high-redshift clusters (z > 0.5) have 
suffered from small number statistics in the simulated sample and the implementation 
of baryonic physics is likely to affect the outcome. We analyse zoomed-in simulations 
of a fairly large sample of cluster-sized objects, with M v ; r > 3 x 10 14 /i _1 M Q , identi- 
fied at z = 0.25 and z = 0.5, for a concordance ACDM cosmology. These simulations 
have been carried out by incrementally increasing the physics considered. We start 
with dark matter only simulations, and then add gas hydrodynamics, with different 
treatments of baryonic processes: nonradiativc cooling, radiative cooling with star for- 
mation and galactic winds powered by supernova explosions, and finally including the 
effect of AGN feedback. Our analysis of strong lensing properties is based on the com- 
puation of the cross-section for the formation of giant arcs and of the Einstein radii. 
We find that the addition of gas in non-radiative simulations does not change the 
strong lensing predictions significantly, but gas cooling and star formation together 
significantly increase the number of expected giant arcs and the Einstein radii, partic- 
ularly for lower rcdshift clusters and lower source redshifts. Further inclusion of AGN 
feedback reduces the predicted strong lensing efficiencies such that the lensing prob- 
ability distributions becomes closer to those obtained for simulations including only 
dark matter. Our results indicate that the inclusion of baryonic physics in simulations 
will not solve the arc-statistics problem at low redshifts, when the physical processes 
included provide a realistic description of cooling in the central regions of galaxy clus- 
ters. As outcomes of our analysis, we encourage the adoption of Einstein radii as a 
robust measure of strong lensing efficiency, and provide the ACDM predictions to be 
used for future comparisons with high-redshift cluster samples. 

Key words: gravitational lensing: strong - galaxies: clusters - methods: N-body 
simulations - cosmology: theory 
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The hierarchical formation paradigm describes a scenario 
in which small structures form earlier and then merge to 
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form more massive objects. This description within the 
ACDM model has been highly successful at explaining 
observations of structure on a larg e range of scales and 
redshifts (e.g. lEfstathiou et al.lll990l ; IWhite fc FrenkHl99ll ; 
iKomatsu et all" 201ll). including the formation of galaxy 
clusters (e.g. iKravtsov fc Borgani|[20ll l . However, the inter- 
nal structure of galaxy clusters has posed a challenge. The 
earliest comparisons between simulated clusters and the ob- 
served frequency of gravitational lensing arcs revealed a se- 
rious disc repancy between the observations and A CDM pre- 
dictions (|Bartelmann et al. I [19981 ; iLi et al]|2005l ). More re- 
cent comparisons to X-ray selected clusters find that the 
discrepancy remains for z < 0.3 and z > 0.5 but results are 
not definitive due to small simulated sample sizes at higher 
redsh ifts and uncertainties wi t h selection procedu res and bi- 
ases (|Meneghetti et al.l [20TTI : iHoresh et all l201lh . In addi- 
tion, observations of clusters at redshifts 0.1 < z < 0.6 have 
found high er concentrations and larger Einstein radii than 
predicted dBroadhurst fc Bark ana 200§j; ISadeh fc Rephaehl 
2008; Zitr in et al. 2011). These are related problems since 
both suggest a failure of ACDM to correctly infer the distri- 
bution of matter in the cores of clusters; the solution may 
be found in the details of the simulations used to formulate 
the predictions. 

Simulations which only incorporate collisionless parti- 
cles are commonplace in the literature. Dark matter struc- 
ture is evolved under gravity with the effect of baryons 
ignored. On large-scales, the distribution of gas follows 
the dark matter potential wells, so the collisionless simu- 
lations will suffice; however on galaxy- and group- scales, 
gas cooling, star formation and feedb ack can alter the 
gravi t ational potent i al sig nificantly (e.g. iBlumenthal et all 
1 19861 ; iGnedin et al.l 120041 ). Early simulations that were 
used to model the cluster lenses d id not incorporate 
sufficient relevant phys ical processes (|Hattori et al.l [l997l ; 
ISadeh fc Rephaelill200Sl} . These simulations were originally 
justified by the premise that the temperature of the ICM 
is too hot to allow efficient cooling, so the baryonic mat- 
ter must simply follow the dark matter potential and 
would not have a significa nt impact on its shape (e.g. 
iBroadhurst fc Barkanal2008l ) . The role of baryons in shaping 
galaxy clusters has recently become a key point of discussion 
since hydrodynamic simulations have become more feasible; 
the findings are that baryons do, in fact, affect the shapes of 
gravitational potential of the simulated haloes. Gas cooling, 
for example, leads to adiabatic compression of galaxy haloes 
preferentially in the cluster centre, allowing more m ass to be 
deposited around the core (|Barkana fc L ocb 2010]). Cooling 
and star formation serve to steepen the density profile of 
clusters and thus increase t he mass concentr ation and the 
strong lensing cross-section dLewis et alj 2000l; 
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2008; Rudd ct al 



200sl ; ICui et al-I^Oll ). Unfortunately, these processes lead 



to the well known overcooling problem, in which the cold gas 
mass and stellar mass in the cluster core is overestimated. 

Including in simulations feedback from gas accretion 
onto supermassive b lack holes (SMBHs) is a promising solu- 
tion for overcooling (ISiiacki et al.ll2007l; pTeyssier et al.l 20 111 ; 
iMcCarthv et al.ll2010l ; iFabian et alj|20ich . while simultane- 
ously reprodjwin^_tlredropin the cosmic SFR for z < 2 
(e.g.. Ivan de Voort et alj|201ll ). The resulting feedback from 
active galactic nuclei (AGN) significantly reduces the gas 



fraction in simulated galaxy-groups and poor clusters (T < 
2keV), but cannot remo ve gas from the deep potential wells 
of rich, massive clusters (IPuchwein et al . 2008; F abian et al.l 
l20ld ; IMcCarthv et alj|2010t ). While the total baryonic con- 
tent of massive clusters is relatively unaffected by cooling, 
star formation and energy feedback, such processes can im- 
ply a quite pronounced redistrib ution of baryons, p articu- 
larly close to the cluster centre. iDuffv et al] (|2010l ) anal- 
ysed simulated haloes at z — ranging from galaxy to 
cluster scales. They demonstrated that w hile radiative coo l- 
ing increases mass concentration(see also lRudd et alj|2008l ). 
AGN feedback has the opposite effect; clusters modelled 
with AGN feedback have mass concentrations equal to those 
mode lled with da r k mat ter only. Consistent with these re- 
sults, iMead et ail (|2010T ) demonstrated that AGN feedback 
is able to reduce the strong lensing cross section for clus- 
ters at z = 0.2. Four of the five clusters they simulated had 
strong lensing efficiencies consistent with their dark-matter 
counterparts, while the fifth remained a stronger lens. 

Many of the previous works investigating the role of 
baryons in modifying cluster cores have told similar stories: 
while cooling and star-formation either directly or indirectly 
leads to an increase of mass in the core, and therefore an in- 
crease in the strong lensing efficiency, AGN feedback negates 
this to some extent. Unfortunately previous studies have 
had limited samples of massive clusters, especially at high- 
redshift z > 0.3, and have utilised simulations generated 
with values for cosmological paramet ers that have since bee n 
revised. Th e stron g len sing studies o f Puchwein et al.l l|2005l ) , 
IRozo et all (|200Sl ) and IMead et all (|2010l ) were undertaken 
in the WMAP-1 typ e cosmology with as = 0.9. As shown by 
IMaccio et al.l (|2008l ), changing the cosmology from WMAP- 
1 to WMAP-3 best fitting models can lead to a 20% decrease 
in the predicted mass concentration of relaxed clusters. This 
is primarily due to the decrease of the power-spectrum nor- 
malisation, as, from 0.9 to 0.7, and the subsequent delay in 
the assembly of haloes. The current WMAP-7 best-fit val- 
ues lie between these two extremes; the expected reduction 
in the concentration of masses considered in t his work from 
WMAP1 to WMAP5/7 is closer to 14% (e.g. IMaccio et ail 
|200S| ; IDuffv et al.ll2008l . but see also lPrada et al.ll2012l )~ 

In this work we improve upon the aforementioned works 
by analysing the strong lensing efficiencies of a larger sam- 
ple of clusters simulated within the currently favoured cos- 
mology (as = 0.8) with only dark matter (DM), and also 
including the hydrodynamical treatment of gas, along with 
cooling, star formation and feedback from both supernova 
(SN) and AGN. Strong lensing properties of our fairly large 
sample of massive clusters are analysed at two redshifts, 
z — 0.25 and 0.5. A description of the cosmological simu- 
lations and the cluster sample follows in Section [2] Earlier 
studies have characterised strong lensing by the cross sec- 
tion for the formation of giant arcs. A comparison with ob- 
servational data requires the cross-section, as a function of 
source redshift, to be convolved with the source redshift dis- 
tribution; uncertainties in this source redshift distribution 
weaken the cosmological test. We, therefore, propose that 
the Einstein radius is a more robust proxy for strong lens- 
ing; critical curves are determined at a single source redshift 
(usually z s = 2), regardless of the range of source redshifts 
used for lens mass reconstructions. Therefore, in Section [3] 
and Section [4] we present the strong-lensing properties for 
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the relaxed sub-sample and discuss the influence of the bary- 
onic processes; in the former section we characterise strong- 
lensing with the cross-section for the formation of giant-arcs, 
while in the latter we consider the Einstein Radii instead. 
At the end of this section, we describe the strong-lensing 
properties of unrelaxed sub-sample. These results are then 
interpreted in terms of variation of the density profiles in- 
duced by the presence of baryons in the simulated clusters. 
Finally, we summarise of our findings in Section [5] 



2 SIMULATIONS 

2.1 The simulated set of clusters 

The set of simulated clusters an alysed in this study have 
been previousl y pres ented by (|Fabian et alj |2010| ) and 
iBonafede et~ai1 (|201lT ). and a more comprehensive descrip- 
tion will be provided in a forthcoming paper (Planelles et 
al. in prep), but we describe them in brief here. Cluster 
halos have been identified in a low-resolution simulation 
box having a periodic co-moving size of 1 h _1 Gpc for 
a flat ACDM model whose cosmological parameters were 
chosen as follows: present day vacuum density parameter, 
f^A.u = 0.76; matter density parameter, Qm,o = 0.24; baryon 
density parameter, Qb,o = 0.04; Hubble constant h = 0.72; 
normalisation of the matter power spectrum as = 0.8; and 
primordial power spectrum P(k) oc k n with n = 0.96. The 
parent simulation followed 1024 3 collision-less particles in 
the box. Clusters were identified at z = using a stan- 
dard Friends-of-Friends (FoF) algorithm, and Lagrangian 
regions around 24 of the clusters found to have masses 
Mfof > lO 18 h -1 M were re-simulated at higher resolu- 
tion employing the Zoomed Initial Conditions code (ZIC; 
iTormen et al.l [1997]) , while resolution is progressively de- 
graded outside these regions, so as to save computational 
time while still providing a correct description of the large- 
scale tidal field. The Lagrangian regions were large enough 
to ensure that only high-resolution particles are present 
within five virial-radii of the central cluster. 

Simulations have been carried out using the TreePM- 
SPH GADGET-3 c ode, a newer version of the original 
GADGET-2 code bv lSpringell (120051 ) that adopted a more ef- 
ficient domain decomposition to improve the work-load bal- 
ance. Each re-simulation has been repeated with a unique 
set of included baryonic processes, some of which are anal- 
ysed here: one simulation-set uses collisionless dark matter 
particles only; one is a non-radiative run; another two imple- 
ment a number of processes associated with baryons; one of 
these two includes AGN feedback. The basic characteristics 
of these re-simulation sets are described here below. 

DM : simulations including only dark matter par- 
ticles, that in the high-resolution region have a mass 
rriDM = 1O 9 /i _1 M0. The Plummer-equivalent co-moving 
softening length for gravitational force in the high- 
resolution region is fixed to epj = 5/i _1 kpc physical at 
2 < 2 while being fixed to ep; = 15fa _1 kpc comoving at 
higher redshift. 

NR. : non-radiative hydrodynamical simulations. Initial 
conditions for these hydrodynamical simulations are gen- 
erated starting from those of the DM-only simulations, and 



splitting each particles in the high resolution region into one 
dark matter and one gas particle, with their masses chosen so 
to reproduce the assumed cosmic baryon fraction. The mass 
of each DM particle is then m D M = 8.47 • 10 s h _1 M and 
the mass of each gas particle is m gas = 1.53 ■ lO 8 h _1 M0. 
For the computation of the hydrodynamical forces we as- 
sume the minimum value attainable by the SPH smoothing 
length of the B-spline interpolating kernel to be half of the 
corresponding value of the gravitational softening length. No 
radiative cooling is involved. 

CSF : hydrodynamical simulations including the effect of 
cooling, star formation and SN feedback. Radiative cool- 
ing rates are computed by f o llowin g the same procedure 
presented by IWiersma et all l|2009l ). We account for the 
the presence of the cosmic microwave background (CMB) 
and for the model of UV/X-ray backgrou nd radiation from 
quasa rs and galaxies, as computed by lHaardt fe Madaul 
(|200ll ). Contributions to cooling from each one of eleven 
elements (H, He, C, C, N, O, Ne, Mg, Si, S, Ca, Fe) have 
been pre-computed usi ng the publicly avail able CLOUDY 
photo-ionisation code (|Ferland et all I1998T ) for an opti- 
cally thin gas in (photo-)ionisation equilibrium. Gas par- 
ticles above a given threshold density are treated as multi- 
phase, so as to provide a sub-resolution description of the 
inter-stell ar medium, according t o the model originally de- 
scribed bv lSpringel fc Hernquistl (|2003l ). Within each multi- 
phase gas particle, a cold and a hot-phase coexist in pres- 
sure equilibrium, with the cold phase providing the reser- 
voir of star formation. We include a detailed description of 
metal production contributed by SN-II, S N-Ia and low and 
interm ediate mass stars, as described by iTornatore et all 
|2007l ). Stars of different mass, distributed according to 
a Chabrier IMF (|Chabrierl [20031 ) ■ release metals over the 
time-scale determined by the corresponding mass-d ependent 
life-times (taken from IPadovani fc Matteuccilll993l ). Kinetic 
feedback contributed by S N-II is implemented ac c ording to 
the scheme introduced by ISpringel fc Hernquistl (|2003l ): a 
multi-phase star particle is assigned a probability to be 
uploaded in galactic outflows, which is proportional to its 
star formation rate. For this set of simulations we assume 
v w = 500 km s^ 1 for the wind velocity. 

AGN : the same as CSF, with a lower wind velocity of 
v w = 350 km s~ , also including the effect of AGN feed- 
back. In the model for AGN feedback, released energy re- 
sults from gas accretion onto SMBHs. This model intro- 
duces som e modifications with r espect to that originally pre- 
sented by ISpringel et all (|2005l ) (SMH hereafter), to which 
is largely inspired, and will be described in detail by Dolag 
et al. (2012, in preparation). BHs are described as sink par- 
ticles, which grow their mass by gas accretion and merging 
with other BHs. Gas accretion proceeds at a Bondi rate, 
while being Eddington-limited. Once the accretion rate is 
computed for each BH particle, a stochastic criterion is used 
to decide which of the surrounding gas particles contribute 
to the accretion. Unlike in SMH, in which a selected gas 
particle contributes to accretion with all its mass, we in- 
cluded the possibility for a gas particle to accrete only with 
a slice of its mass, which corresponds to 1/4 of its original 
mass. In this way, each gas particle can contribute with up 
to four "generations" of BH accretion events, thus providing 
a m ore continuous desc ription of the accretion process (see 
also iFabian et afll201Ch . BH particles are initially seeded 
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with a mass of 0.05 moM.io, where moM.io is the DM par- 
ticle mass in units of lO lo h _1 M0. Seeding of BH particles 
takes place in halos when they first reach a minimum friend- 
of-friend (FoF) mass of 2.5 x 10 3 hidm.io (using a linking 
length of 0.16 in units of the mean interparticle separation in 
the high-resolution region), with the further condition that 
such halos should contain a minimum mass fraction in stars 
of 0.02. The first condition on the minimum halo mass guar- 
antees that such halos are resolved with at least ^ 200 DM 
particles, while the second condition requires that substan- 
tial star formation took place in such halos. This criterion 
prevents seeding BHs in halos possibly located at the border 
of the high resolution region, which spuriously contain a low 
amount of cooled gas, due to the interaction with nearby 
low-resolution DM particles. Eddington-limited Bondi ac- 
cretion produces a radiated energy which correspond to a 
fraction e r = 0.1 of the rest-mass energy of the accreted gas, 
which is determined by the radiation efficiency parameter t r . 
The BH mass is correspondingly decreased by this amount. 
A fraction of this radiated energy is thermally coupled to 
the surrounding gas. We use e/ = 0.05 for this feedback 
efficiency, which increases to e/ = 0.2 whenever accretion 
enters in the quiescent "radio" mode and takes place at a 
rate smaller than one-hundredth of the Eddington limit (e.g. 
ISiiacki et al1l2007l ; lFabian et alj|20ld ). 



2.2 The identification of clusters 

The cluster haloes are identified as follows. Firstly, a stan- 
dard Friends-of-Friends algorithm is run over the DM par- 
ticles in the high-resolution regions, using a linking length 
of 0.16 in units of the mean- interparticle separation. Within 
each FoF group, we identify the position of the particle with 
the minimum gravitational potential, which is then taken as 
the centre from where clusters are then identified accord- 
ing to a spherical overdensity (SO) method. The virial ra- 
dius is defined as the smallest radius of a sphere centred 
on the cluster, for which the mean density falls below the 
virial overdensity. The virial overdensity is measured relative 
to the critical density and calcu lated using the fitting for- 
mula of lBrvan fc Normanl (|l99ct ). so for clusters at z = 0.25, 
A c £3 112 and for clusters at z — 0.5, A c ~ 129. 

Cluster haloes are chosen for our strong lensing anal- 
ysis only if no low-resolution particles contaminate the re- 
gion within five virial radii of the cluster centre. In Figure [T] 
we show the number of clusters above a given mass limit 
at the two redshifts considered in our analysis. Note that 
this is not equivalent to a halo mass function since the re- 
simulated regions do not cover the entire box. The sample 
has 42 clusters with M vir > 3 x 10 14 h~ 1 M Q at a redshift 
of z = 0.25. Moving to higher redshifts, we find 34 clusters 
with Mvir > 3 x 10 14 ft -1 M Q at a redshift of z = 0.5. 



2.3 Relaxed clusters 

If mergers between clusters take place along the line of sight 
the strong lensing efficiency can be enhanced. If the merger 
occurs across the sky, we will see multiple critical curves, 
which positively biases the cross section for giant arcs. In 
the case that critical curves are merging, both the Einstein 
radii and the cross section for giant arcs would be biased 
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Figure 1. The cumulative virial-mass distribution of clusters in 
the sample at redshifts z = 0.25 (orange dotted curve) and z = 0.5 
(purple dashed curve). 



large. This could confuse our interpretation of the effects of 
baryon physics. 

In order to characterise the degree of relaxation of each 
cluster, we compute the quantity s max = {max(sQ) : 0.05 < 
£ < 2}, where 



|»'COm(< CR-vir) — rMinPot| 
CR-vir 



(1) 



In the above equation, rMinPot is the position of the particle 
with the minimum gravitational potential and rcoM is the 
centre of mass of all the matter within £R v ir. In this way, S( 
represents the offset between these two measures of "centre" 
for any aperture radii, CR-vir, so that s max is the maximum 
offset over a range of aperture radii, which are found by 
allowing £ to vary from 0.05 to 2 in 30 logarithmic steps. We 
define a relaxed cluster as that for which the counterpart 
in the DM simulation has a maximum offset parameter of 
s < 0.1. 

As a visual aid, the offset parameter for each cluster 
at z = 0.5 is shown in Figure [2] plotted against the aper- 
ture radius; each curve corresponds to an individual clus- 
ter. Any cluster for which the curve exceeds the dotted line 
(sf = 0.1) for any radius in the given range is considered 
to be unrelaxed. Our adopted criterion to classify relaxed 
and unrelaxed clusters corresponds to similar definitions in 
the literature when we choose only one a perture radius, the 
virial radius i.e. when we set C = 1 (e g.lCrone et al.lll996l; 
Thomas et ai1ll998l; iNeto et alj|2007t iD'Onghia fc Navarro! 
20071 ; I Power et alj|201lh . However, we note that adopting 
this common choice we would be less sensitive to complex 
mass distributions near the core of a cluster, which is more 
important in the context of strong lensing. Our definition, 
in a sense, weights the disturbance caused by substructure 
by the inverse of radial position of the substructure. In order 
to verify whether there is any correlation within our sam- 
ple between degree of relaxation and cluster mass, we plot 
in Figure [3] the positions of our clusters in the s max -M v i r 
plane. There is negligible correlation between the cluster's 
mass and the offset parameter that describes the degree of 
relaxation. The choice of threshold for s max above which a 
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Figure 2. The offset parameter, sr, as a function of aperture 
radius within which centre-of-mass is calculated, for all clusters 
at z = 0.5. The dotted horizontal line shows the value of the 
maximum offset parameter, s ma x = 0.1, allowed for a cluster to be 
classified as relaxed. Red and blue curves correspond to clusters 
which are then classified as unrelaxed and relaxed, respectively. 
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cluster is considered unrelaxed does not appear to bias the 
sample in terms of mass. The total sample of 42 clusters at 
2 = 0.25 includes 17 relaxed clusters, and of the 34 clusters 
in total at z = 0.5, 14 are classified as relaxed. In Section [3] 
and at the beginning of Section [4] we present the statistical 
results only for the relaxed sub-sample. In Section 14.31 we 
compare these results to those for unrelaxed clusters. For 
comparison, we note that if we had defined the relaxed sub- 
sample as those for which the offset at the virial radius is 
s\ < 0.07, as is common in the literature, we would have 
23 relaxed clusters at z = 0.25 and 12 relaxed clusters at 
z = 0.5. 



3 CROSS SECTION FOR GIANT ARCS 

The cross section for giant arcs is defined as the area in the 
source plane in which if a distant galaxy (at that redshift) 
was located, it would appear as a highly elongated giant 
arc due to the effects of strong gravitational lensing. Every 
line of sight through each simulated cluster will, in general, 
provide a unique cross section. We now describe the proce- 
dure used to calculate the giant arc cross section for each 
line of sight analysed for each cluster in our simulated sam- 
ple. Throughout the present work, we refer to gravitational 
lensin g quantities following the notation of ISchneider et all 
(|l992t ). The lensing mass is chosen to be contributed by all 
particles within two virial radii of the cluster centre. We 
apply the thin lens approximation: the assumption that the 
lensing mass is constrained to a single plane. This treatment 
is valid for the present study, since the sizes of the cluster 
lenses are much smaller than D s , D d and Dd s - the angu- 
lar diameter distances from the observer to the source, from 
the observer to the lens, and from the lens to the source, 
respectively. The convergence, k, can be defined by: 

«=^, (2) 

where E is the surface density and the critical surface density 
for gravitational lensing is given by 

_ c 2 D s 
^ - 4nG D d D d / [6) 

Note that from here onwards, the redshift of the cluster 
lenses is denoted zl, while the redshift of the background 
source galaxies is denoted z s . 



3.1 Calculating the lensing cross section 

In order to calculate the lensing cross-section, we first mea- 
sure the deflection of light-rays from background source 
galaxies across a field of view in which we expect to search 
for giant arcs. Each ray of light is deflected at the lens plane; 
the deflection angle is related to the convergence by: 



a(x) 



d x'k(x.') — 
l x 



(4) 



where x is a dimensionless position in the lens plane. It is 
possible to choose the scaling such that the gravitational 
lens equation is given by: 



13 = 9- a(0), 



(5) 



where (3 is the angular source position and is the angu- 
lar position of the image on the sky. Producing maps of the 
deflection angle at each grid point requires a few steps, as 
shown in the cartoon diagram in Figure [4] We describe the 
procedure below. The angular and spatial resolutions are 
quoted for sources at z s = 2 and clusters at zl = 0.25, 
but quantities for clusters at zl = 0.5 are given in square 
brackets. The lensing mass is projected onto one of two pos- 
sible lens planes, depending on the position of the particles. 
Those that lie at a projected co-moving distance of more 
than 0.9 h _1 Mpc [1.35 h _1 Mpc] from the cluster centre are 
placed on a low-resolution 2048 x 2048 grid. The angular res- 
olution of this grid is 1.7 arcsec [0.9 arcsec]. Particles that 
are projected closer to the cluster core, and are therefore 
responsible for the bulk of the strong lensing, are placed on 
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Figure 4. A cartoon diagram describing the calculation of deflection angle, a, across a high resolution grid. The lensing mass is divided 
into two parts: an inner and an outer region (shown in the bottom and top rows, respectively), each of which are projected onto a 
small high-resolution and large low-resolution grid, respectively. The surface density distribution on each grid, k, is smoothed, then the 
deflection angle is calculated using equation [4] The total deflection angle is the sum of these two. However the deflection angle calculated 
in the inner region of the large grid must first be interpolated onto the small high-resolution grid. See text for further details 



a high-resolution 2048 x 2048 grid. The angular resolution 
of this grid is 0.3 arcsec [0.2 arcsec]. The mass on each grid 
is smoothed with a truncated Gaussian filter of size a — 5 
h _1 kpc, to match the force-softening length of the simula- 
tion. The size of each grid is sufficiently larger than the pro- 
jected extent of the mass on that grid, so that smoothing 
does not result in any loss of mass. 

We convolve each of the two convergence maps with the 
appropriate kernel from equation [4] This leaves us with two 
maps of the deflection angle: one at high-resolution where 
the deflection is solely due to the mass in the projected clus- 
ter core; the other at low-resolution from the rest of the lens- 
ing mass. For the subsequent ray-shooting procedure, we use 
a single map of deflection angle for each cluster projection; 
the final map has the same size and resolution as the high- 
resolution grid described above. The deflection angle at each 
grid point is determined with bilinear interpolation from the 
closest grid-points of the low-resolution map, and adding the 
corresponding grid-point from the high-resolution map. The 
resulting deflection angle map has an angular resolution of 
0'.'3 for clusters at zl = 0.25 and 0"2 for clusters at zl = 0.5. 

Calculating the cross-section for giant arc formation re- 
quires a ray-shooting and arc-ide ntifying procedure , which 
has been previously p r esente d in Mcncghc tti et all (|2000l ) 
and iMeneghetti et aLl (120051 ), Elliptical sources with an 
equivalent radius of size 0"5 are placed throughout the 
source plane; more sources are placed in regions identified 
as caustics. The deflection angle maps allow mock images 
to be generated from the sources. The images are fitted 
to ellipses and the length-to-width ratio, L/W, is deter- 
mined. If L/W surpasses some threshold elongation, rj, the 
image is identified as a 'giant arc'. The cross-section for gi- 
ant arc formation, a v , is the area in the source plane cov- 
ered by sources that are mapped onto giant arcs. For the 
present study, the len gth-to-width thr e shold i s chosen to 
be r\ = 7.5 — following iPuchwein et all (120051 ). iMead et all 
(|2010l ) and IMeneghetti et al.l ( 201ll ) — because much larger 
cross-sections are subject to small number statistics, and 



much smaller cross-sections are too sensitive to the intrinsic 
cllipticities of the source galaxy. We discuss other choices of 
threshold of Section \3. 2 1 

We present results in this section and Section [3] for 
sources at redshift z s = 2, but include the statistical results 
for a source redshift of z s = 1 in Table 1. For these lower 
source redshifts, we would expect smaller Einstein radii, so 
the deflection angle maps have a higher angular resolution: 
0"2 for clusters at zl = 0.25 and O'.T for clusters at zl = 0.5. 



3.2 The impact of baryons on the cross section 
for giant arcs 

The baryonic processes as implemented in the hydrody- 
namical simulations affect the details of the structure of 
each cluster halo. Figure [5] shows a single cluster of mass 
M V ir = 8 x 10 14 h~ 1 M e in our zl = 0.5 sample as found 
in each of the four re-simulations. The basic structure is 
recognisable in each simulation and the subhaloes appear in 
much the same positions, as seen in the convergence maps 
along the top row. However, the distribution of matter near 
the clusters centre is sufficiently altered by baryonic physics 
so that the caustics, as shown in the bottom row of the 
same figure, are distinct. In particular, note that the CSF 
simulation produces larger source caustics, shown in blue, 
and larger tangential image caustics — or critical curves 
— shown in red. The tangential part of the image caustic 
extends around a substructure that was sub-critical in the 
other simulations. For this particular cluster, as seen from 
one line of sight, the strong lensing efficiency is highest when 
the cluster is generated within the CSF simulation. 

It is our aim to determine whether this is generally true 
for simulated galaxy clusters, and in particular, how the 
strong lensing properties predicted by the AGN simulations 
compare with the other simulations. We have a large sample 
of galaxy cluster, but we should not consider only one line of 
sight through each cluster. Due to the aspherical nature of 
clusters, as well as the presence of substructure, each line of 



Baryons & Cluster Strong Lensing 7 





arcsecs arcsecs arcsecs arcsecs 

Figure 5. One cluster at zj, = 0.5 seen from the same viewpoint within different simulations: DM (left panel), NR (second panel), CSF 
(third panel), and AGN (right panel). Along the top row, we show the convergence map of the cluster, ~ 410 arcsecs across, in which 
brighter region represent higher surface densities; along the bottom row is the corresponding source caustic in blue and image caustic, 
a.k.a the critical curve, in red, both determined for a source at redshift z s = 2 




Figure 7. The cumulative probability distribution for the source-plane cross-section for the formation of giant arcs combining 50 lines 
of sight through the relaxed clusters at = 0.25 (left panel) and z^ = 0.5 (right panel). Different lines and colours have the same 
meaning as in Figure [6] The arrows mark the median values of the 1x7.6 cross section for the models. 



sight produces a unique cross section. Therefore, each clus- 
ter is associated with a distri bution of possible cross sections 
(see, for example, figure 3 of lDalal et al.|[2004h . In Figure [U 
we show such a distribution by measuring the cross section, 
(T7.5, for 50 lines of sight through one cluster in our sample, 
also seen in Figure [5] In the present study, we analyse pro- 
jections of each cluster along 50 randomly chosen lines of 
sight. The results for all our sets of simulated clusters are 
combined in the left-hand panel of Figure [7] 

We conduct a two-sample Kolmogorov-Smirnov test to 
compare the probability distribution of (77.5 values for each 
simulation against the same distribution for the DM simula- 



tions. The D-statistic is the maximum difference between the 
two cumulative probability distributions. Ideally we would 
like to determine the probability, p a , that the D-statistic 
would be atleast as large as that measured assuming the 
samples are drawn from the same underlying distribution. 
We recognise that while the clusters are independent, the 
many lines of sight analysed for each cluster are not inde- 
pendent from each other. In an actual observational sample, 
one would measure a single strong-lensing efficiency for each 
cluster; however using a simulated sample, scatter due to 
orientation has to be taken into account. Therefore, stack- 
ing the results for all lines of sight analysed, we calculate 
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Figure 6. The probability distributions for the cross section that 
would be measured for a single simulated cluster at Zi, = 0.5 for 
a source redshift of z s =2; the cluster is the same as that shown 
in Figure [5] Different curves refer to the distributions produced 
by the four different simulations: DM (black solid line), NR (blue 
dashed line), CSF (green dash-dotted line). AGN (red dotted line) 



the D-statistic. Using the total number of lines of sight (50 
times the number of clusters) as the 'sample size' provides 
a lower limit for the p- value, p™ m ■ Alternatively, by letting 
the number of clusters be the 'sample size', we can obtain 
an upper limit on the p- value, p™ ax . Both are listed in Table 
1. If the clusters were all perfectly spherically symmetrical, 
so that all lines of sight resulted in the same value of 0-7.5, 
then p™ ax would be equivalent to the required p a . 

Table 1 also include the typical cross-sections for clus- 
ters in each simulation, for low and high source redshifts: 
z 3 — 1 and z 3 — 2. These are calculated by taking the me- 
dian value of the cross-section over the 50 lines of sight anal- 
ysed for each clusters, then averaging over all clusters. Com- 
paring the results for DM and NR simulations, we find that the 
ability to lens background galaxies into giant arcs is similar, 
independently of whether the cluster lens is simulated with 
collision-less particles only, or with additional non-radiative 
gas. As seen for the CSF simulation, when cooling and star- 
formation are implemented, the cross-section for the forma- 
tion of giant arcs is boosted by a factor of ~ 1.5 for high 
source redshifts (z s = 2) and ~ 2 for lower source redshifts 
(z 3 — 1). This is a slightly small er 'boost' rel a tive t o the 
findings of IPuchwein et all (|2005l ). iRozo et all (|2008l ) and 
iMead et all (|2010l ). This difference with respect to previ- 



ous simulation models can be understood in terms of the 
more efficient SN feedback implemented in our CSF sim- 
ulation set, which more efficiently countera cts the effect of 
halo c ontraction induced by cooling. In fact. lPuchwein et al.l 
|2005l ) use the same model for galactic ejecta, but with a 
smaller veloc ity of galactic win ds, with v w = 350 km s _1 . 
Furthermore. IRozo et all (|2008l ) and lMead et all (|2010D only 
include thermal schemes for SN feedback, which are rather 
inefficient in regulating cooling at the centre of cluster-sized 
halos. 

AGN feedback reduces the giant-arc cross-section and 
makes it comparable to the predicted cross-section from 
collision-less simulations. When we compare clusters mod- 
elled with dark matter only (DM results shown in black) and 
their counterparts modelled with the complete set of bary- 



onic processes, including AGN feedback (AGN shown in red), 
we find surprisingly little difference in the typical cross- 
sections for clusters at zl = 0.25 - now a boost of only 
~ 10 per cent. However, at zl = 0.5 clusters modelled 
with AGN feedback are ~ 30 per cent more efficient at 
lensing background galaxies into giant arcs. Given the re- 
sults for Einstein radii (see Section T4.2I below) . this is likely 
due to the presence of substructures that produce additional 
arcs rather than the re-distribution of mass near the cluster 
centre. The difference in the results between CSF and the 
AGN simulations arise from a combination of the reduced SN 
wind-speed and the introduction of thermal AGN feedback. 
However, if the CSF simulations had a similarly reduced 
wind-speed, the reduced feedback efficiency would lead to 
increased cooling and a higher lensing cross-section, there- 
fore it is the AGN feedb ack that is r espon sible for reduc- 
ing the lensing efficiency. IMead et al.l (|2010l ) analysed three 
projections for each of five clusters at zl = 0.2 with simi- 
lar simulation sets. For four of these five clusters the strong 
lensing cross-sections for the dark matter simulation and 
their AGN simulation were very similar. At this redshift, 
given the small sample and number of projections, this is 
consistent with our findings. 

As described earlier, images are identified as arcs if their 
length-to-width ratio surpasses the threshold elongation, r\. 
Although extreme values for this threshold are not useful, 
the exact value is somewhat arbitrarily chosen; to be sure 
that the results are not qualitatively affected by this choice, 
we plot the cross section as a function of r\ in Figure [5] The 
lines represent the median cross section over the 50 lines of 
sight for each individual cluster, then averaged over all the 
clusters in the (relaxed) sample. The error bars represent 
the 16th and 84th percentile for individual clusters, then 
averaged over all the clusters. These error-bars, therefore, 
reflect the triaxiality of the clusters. Unsurprisingly, a larger 
and more extreme choice of threshold allows a smaller cross 
section, since less sources would be able to produce such 
long and thin arcs, but the qualitative differences between 
the different simulation sets are the same. 

By comparing the error -bars in Figure [§ ] for all re- 
simulations, we find that as in lRozo et al.l (|2008l ). the spread 
in the distribution of possible cross-sections is also reduced 
significantly in the CSF simulations. These clusters are less 
susceptible to light-of-sight effects than those modelled with 
dark matter only or non-radiative gas, predominantly due 
to the more spherical shape of the clusters resulting from 
the condensing of baryons and t he subsequent response 
from the dark matt e r component dBlumenthal et alj|l98a: 
Gnedin et al. I |2004 iKazantzidis et alj |2004| ; iBrvan et"ail 



20121 ). The steeper inner profile (see Section \^A\ makes the 



clusters less susceptible to coincidental mass along the line of 
sight. Although AGN feedback regulates the isotropic con- 
densation of baryons, the reduced spread in cross-sections is 
also found for the AGN clusters, particularly at zl = 0.5. 

Figure [9] shows the giant arc cross-section, 07. 5, plot- 
ted against the mass of each clusters. Each cluster is repre- 
sented as a single point in the graph, with different colours 
representing the counterparts in the four re-simulations. 
The error-bars mark the 16th and 84th percentiles of 07. 5 
found over the 50 lines of sight analysed. The figure sug- 
gests that there exists a scaling relation between lensing 
efficiency and mass, but that the correlation is tighter at 
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Figure 8. The giant arc cross section plotted against chosen elongation threshold r\ for the relaxed cluster sub-sample at zl = 0.25 (left 
panels) and zl = 0.5 (right panels). The cross section for a given value of r\ is the median of 50 lines of sight through each cluster, then 
averaged over all clusters. The errorbars mark the 16th and 84th percentile for each cluster, then averaged over all clusters. The source 
redshift is z s = 2 
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Figure 9. The cross-section for the formation of giant arcs at source redshift of z s = 2 as a function of cluster mass, for our relaxed 
cluster sub-sample. Each dot represents the median cross-section from the 50 lines of sight analysed for a single cluster, while the error 
bars mark the 16th and 84th percentiles. In the left-hand panels, the characteristic mass is Afsoo while for the right-hand panels, the 
characteristic mass is A/2500 • The clusters at zl = 0.25 are shown on the top row while the clusters at zl = 0.5 are shown on the bottom 
row. 
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Table 1. Statistical Analysis of Strong Lensing Properties of Relaxed Clusters. Column 1: simulation set. 
Columns 2 and 3: redshifts of the lens, zl, and of the source, z s . Column 4: Pearson correlation coefficient 
r for the log(<T7.5)-log(0E) relation. Column 5 and 6: least-squares fit to the log(o-7.5)-log(#E) relation (see 
Eqn. |6j- Column 7: typical value of the Einstein radius, 8^, computed as the median value over 50 lines of 
sight through each cluster, then averaged over all clusters (units of arcseconds). Columns 8 and 9: values 
of p™ n and p™ ax , defined as the lower and upper limits, respectively, of a true p- value from a KS-test 
performed to compare the $e probability distribution of each simulation set against the corresponding result 
for the DM simulation. Column 10: the typical cross section for the lensing of galaxies at z s into giant arcs 
with elongation r\ = 7.5, 177.5, defined as the median value over 50 lines of sight through each cluster, then 
averaged over all clusters (units of 10~ 4 /i~ 2 Mpc 2 ). Columns 11 and 12: the same as in columns 8 and 9, 
respectively, but to compare the 0-7 5 probability distribution. 



Simulation 




Zs 


r 


a 


b 


8e 


Pe 




07.5 


if a 




DM 
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NR 
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< 10- 5 
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1.00 
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20 






8.2 






NR 


0.5 


2 
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0.5 
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4 x 10~ 4 


1.00 
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10 
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higher-overdensities; this reflects the region most responsi- 
ble for strong lensing. The higher scatter when 07.5 is plot- 
ted as a function of M500 reflects the fact that different re- 
simulations create differing distributions of baryons in the 
cluster core, when comparing clusters at a fixed mass at 
lower overdensities. 



4 EINSTEIN RADII 

The cosmological test based on arc statistics is subject to un- 
certainty in the characteristics of the source population. Per- 
forming a comparison with observed arc statistics requires 
one to convolve the predicted giant arc cross-section with 
an assumed redshift distribution for the background galax- 
ies. The uncertainties in this redshift distribution creates an 
additional uncertainty in the cosmological test. This addi- 
tional uncertainty is unnecessary since one might instead 
characterise strong lensing efficiency by the angular scale 
which separates highly magnified images; this is known as 
the Einstein radius. Measuring the Einstein radii requires 
the observer to: measure the shape and redshift of a rea- 
sonable number of highly magnified high redshift sources, 
reconstruct the lens mass distribution, and infer the critical 
curves at a single redshift. This is achievable even if the var- 
ious sources are at different redshifts, and bypasses the need 
to independently measure the expected number density and 
redshift evolution of the sources. We therefore measure the 
Einstein radii for our cluster sample and compare the results 
between the different re-simulations. 



4.1 Measuring the Einstein Radius 

The angular separation of highly magnified background 
galaxies has a formal definition, which is strictly applicable 
only in the case of axially symmetric lenses. However, galaxy 
clusters are not axially symmetric in general; there exist a 
number of different definitions in the literature for how one 
might calculate the Einstein radius for more realistic lenses. 
For example, one might divide the area enclosed within the 
tangential critical curve by n, as per the so-called 'equivalent 
Einstein radius' definition used , e.g.. bvlPuchwein fc Hilbertl 



|2009l ). lzitrin et al l J201ll) andlRedlich et all |20T2T ). On the 
other hand, iBroadhurst fc Barkanal (|200&t ) define the 'effec- 
tive Einstein radius' as the radius that encloses a mean sur- 
face density equal to the critical surface density n — 1 (see 
iPuchwein fc Hilbertl2009l for further discussion on how these 
two d efinitions compare ) . We , instead, follow the defini- 
tion of iMeneghetti et al.l (|201 lh and characterise our statis- 
tics by a 'median Einstein radius': the median distance 
of the tangential critic al points from the clusters centre. 
IMeneghetti et all (|201ll ) has demonstrated that the median 
Einstein radius has a tighter correlation with strong lensing 
cross section than the 'equiva lent Einstein radius'. W hen 
measuring the Einstein radius, IMeneghetti et alj (|201lh de- 
fine the cluster centre to be the location of the maximum 
of the projected mass distribution. In the present work, we 
define the centre to be the projected position of the particle 
in the simulated cluster with the lowest potential. This is 
done in order to avoid attributing the position of the centre 
to a subhalo in cases in which such a subhalo is fortuitously 
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projected so it has a higher projected mass density than the 
centre of the main halo. 

We use the high resolution deflect ion- angle map (as de- 
scribed in Section [3.f I) to identify tangential critical points 
within the same field of view at the same angular resolu- 
tion. There can be complications in the presence of sub- 
structure, which require one to discard some of the critical 
points before measuring the Einstein radius. We show an un- 
relaxed cluster in Figure [TO] as an example. There is a large 
substructure present and, from some lines of sight, visible 
within the field of view. For each line of sight analyzed, crit- 
ical points are mapped out across the lens plane forming a 
'critical curve' also known as an image caustic; these corre- 
spond to the caustics in red. Then, critical points associated 
with the tangentially sheared images (as opposed to radially 
sheared) at the 'primary' critical curve around the centre of 
the cluster are identified, so that 'secondary' critical curves 
associated with substructure do not bias the measurement. 
These are the 'cleaned' curves shown in black. The projected 
distance of each 'clean' point from the cluster centre is mea- 
sured; the Einstein radius, #e, is the median value over all 
these points. 

If the secondary critical curve merges with the primary, 
as shown in the bottom row of Figure llOl they cannot be dis- 
tinguished, and 9e is artificially increased. The cross-section 
""7.5 gives, to a first approximation, the region enclosed by 
the source caustics (shown in blue). The measurement of (T7.5 
can be artificially enhanced by the presence of substructure. 
We discuss this further in Section [473] For a relaxed cluster, 
there would be fewer lines of sight in which the secondary 
and primary critical curves merge, and the cluster will have 
smaller secondary critical curves which have a less significant 
impact on the measurement of (J7.5 and (9e- 

Figure [TT] shows the giant-arc cross section and the Ein- 
stein radi i for corresponding lines of sight through each clus- 
ter. As in lMeneghetti et al.l |201ll ). it is clear that the corre- 
lation is very strong for our relaxed cluster subsample. We 
perform a least-squares fitting to a function of the form: 

log(<77.5) = alog(6> E ) + b, (6) 

and measure the Pearson correlation coefficient for each sim- 
ulation separately, including only lines of sight with non-zero 
values for both the cross section and Einstein radii. The re- 
sults for a, b and r are summarised in Table I. The correla- 
tion coefficient, r, ranges from 0.92 to 0.97 for the relaxed 
clusters at zl = 0.25, and from 0.95 to 0.96 for the relaxed 
clusters at zl = 0.5. Two of the clusters in the NR simula- 
tion are responsible for highly elliptical critical curves and 
elongated source caustics, which results in the scatter above 
the line of best-fit in the left panel of Figure 1111 and makes 
the line of best-fit shallower; this is indicative of the sensi- 
tivity of the fit to su bstructure. Similarly, a recent study by 
iRedlich et all (|2012l ) has concluded that the presence of sub - 
structure and cluster mergers in the lMeneghetti et alj (|201ll ) 
simulated clusters sample — as opposed to semi-analytic 
smooth triaxial cluster-halo models — results in shallower 
slope and higher normalisation. The line of best-fit to Eqn.[S] 
for the NR simulation at zl = 0.5 and z 3 = 2 is the most ap- 
propr iate for comparison to the results of iMeneghetti et al.l 
(|201ll ). which are obtained using the z > 0.5 clusters found 
in the non-radiat i ve Mar eNostrum simulation. Compared to 
IMeneghetti et all (|201lf ) we find a lower normalisation, a, 



and a steeper slope, 6, of the log(<77.5)-log(#E) relation. This 
can be attributed to our deliberate exclusion of clusters with 
significant substructure and merger-activity in our relaxed 
subsample. 

We can expect the slope of the log(a7.5)-l og(^E) relation 
to re flect the inner mass profile of the lens (|Redlich et al.l 
120121) . Dark-matter h aloes, such as those in DM simulations, 
follow NFW profiles (|Navarro et al.lll996l ). with inner den- 
sity profiles that fall off as p(r) oc r _1 . In clusters with 
significant gas cooling, such as those in the CSF simulations, 
the density profile is close to isothermal (p oc r~ 2 ), and 
so, produce a large tangential magnification relative to ra- 
dial magnification, and therefore a boost in 07.5 compared 
to the DM clusters. Since the positions of the critical points 
are relatively stable to the inner slope — assuming no sig- 
nificant redistribution of mass to larger radii — there is a 
greater boost in the giant arcs cross section as opposed to 
Einstein radius, part icularly for low mass lenses (see fig- 
ure 5.7 in lQgurill2004T l. which results in a shallower slope in 
the log(<T7.5)-log(6 | E) relation, and a smaller value of a. For 
low source redshifts, lensing efficiencies are lower, resulting 
in greater scatter as the source sizes approach the resolution 
of the lensing maps; the line of best-fit is sensitive to this 
scatter and less sensitive to the inner mass profiles. 

4.2 The impact of baryons on the Einstein radii 

We plot the cumulative probability distribution for Einstein 
radii in Figure [12] combining all 50 lines of sight through 
each of the relaxed clusters. As in Section |3. 2 1 we conduct a 
two-sample Kolmogorov-Smirnov test to compare the prob- 
ability distribution for each simulation against the DM sim- 
ulation. We want to determine the probability, pg, that the 
D-statistic would be atleast as large as that measured as- 
suming the samples are drawn from the same underlying 
distribution. As noted before, while the clusters are inde- 
pendent, the many lines of sight analysed for each cluster 
are not independent from each other. However, for the sim- 
ulated sample, scatter due to orientation has to be taken 
into account. Stacking the results for all lines of sight anal- 
ysed, we calculate the D-statistic. Defining the 'sample size' 
as the total number of lines of sight (50 times the number 
of clusters) provides a lower limit for the p- value, pg un . Al- 
ternatively, by letting the number of clusters be the 'sample 
size', we can obtain an upper limit on the p- value, p^ ax . 
Both are listed in Table 1. 

Table 1 also includes the typical Einstein radii for clus- 
ters in each simulation, for low and high source redshifts: 
z a — 1 and z s = 2. These are calculated by taking the 
median value of the Einstein radius over the 50 lines of 
sight analysed for each clusters, then averaging over all clus- 
ters. The strong lensing efficiency of a cluster is sensitive 
to the enclosed mass within #e- Therefore, lens-source con- 
figurations that produce larger Einstein radii are less sus- 
ceptible to the effects of baryons, since AGN lead to the 
re-distribution of mass within reasonable small radii; we 
leave a more detailed discussion to Section l4.4l . Comparing 
the results for DM and NR simulations, we find that Einstein 
radii are similar whether the cluster lens is simulated with 
collision-less particles only, or with additional non-radiative 
gas. Comparing the results for DM and CSF simulation, gas 
cooling and star formation increase the predicted Einstein 
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Figure 10. The surface density map (left-hand panels) and image and source caustic (right-hand panels) of an unrelaxed = 0.5 cluster 
in our sample. In the left-hand panels, brighter regions correspond to higher surface density and the image caustic is overlaid in red. In 
the right-hand panels, the image caustic is shown in red, with the 'cleaned' caustic - used to infer the Einstein radius - drawn in black, 
and the source caustic at z s = 2 shown in blue. The 'primary' image caustic is associated with the centre of the cluster, while a large 
substructure present within the field of view produces a 'secondary' image caustic. Seen from one line of sight (top row), the secondary 
caustic is distinct, so the cleaned caustic captures only the primary caustic, while from the other line of sight (bottom row) the primary 
and secondary caustics are merged, so the cleaned caustic captures both, and artificially enhances the inferred Einstein radius. 
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Figure 11. Relationship between the giant-arc cross section, (77,5, and the Einstein radius, 6>e, for each of the 50 orientations of each 
cluster. On the left panel, we show the results from the clusters at zl = 0.25 and on the right panel, we show the results from the clusters 
at = 0.5; both are for source redshift of z s = 2. Results for the four simulations are combined: dark matter only (black); dark matter 
and non-radiative gas (blue); with cooling and star- formation (green); and with AGN feedback (red) 



radii by ~ 10 per cent for high source redshifts (z 3 = 2) and 
~ 20-40 per cent for lower source redshifts (z s = 1). This 
effect is more significant for lower source redshifts, for which 
the Einstein radii are typically smaller. With the inclusion of 
AGN feedback, the typical Einstein radii are reduced with 
respect to the CSF case. For Einstein radii calculated for 
source redshifts of z s = 2, the fractional increase relative 



to the clusters in the DM simulations is only ~ 5 per cent. 
This is the source redshift that reconstructed critical curves 
are commonly scaled to in observational studies. The great- 
est difference between predicted Einstein radii from DM and 
AGN simulations occurs for high-z clusters and low-z sources, 
which corresponds to the smaller Einstein radii (9 < 20"). 

Figure [13] shows the Einstein radius for each cluster as 
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Figure 12. The cumulative distribution of Einstein radii obtained by combining 50 lines of sight through each of the relaxed clusters 
in the sample. The source redshift is assumed here to be z s = 2. On the left panel, we show the results for the 24 relaxed clusters with 
M vir > 3 X 10 14 h~ 1 M e at redshift z L = 0.25. On the right panel, we show the results for the 18 relaxed clusters above the same mass 
limit at redshift z^ =0.5 
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Figure 13. The Einstein radius for a source redshift of z B = 2 versus cluster mass combining 50 lines of sight through the relaxed 
clusters at z^ = 0.5. For the left-hand panels, the characteristic mass is M500 while for the right-hand panels, the characteristic mass is 
M2500 



a function of its mass. We plot the median value of 8e mea- 
sured over the 50 lines of sight for an individual cluster; the 
error bars, denoting the 16th and 84th percentiles, reflect 
a measure of spread associated to the line-of-sight variance. 
We note that one of the clusters has a unusually large scatter 
in Be associated with the different orientations. This clus- 
ter has a small Einstein radius and is probably on a similar 
scale to angular resolution. The cluster mass is measured at 
two different overdensities: A = 500 and A = 2500. The Ein- 
stein radius increases with cluster mass, as expected, but the 
correlation is tighter at higher overdensities. This is in line 
with the expectation that high overdensities are responsible 
for strong lensing. The Einstein radius for a fixed value of 
M500 is slightly larger for clusters in the CSF simulations due 
to the higher concentration of these clusters. The aspheric- 
ity, reflected in the size of the error bars, is also smaller for 
clusters in these simulations. 



4.3 Unrelaxed Clusters 

Unrelaxed clusters produce highly complex caustic struc- 
tures which are the consequence of high mass subhaloes 
that lie within the field of view. These can create non-trivial 
complications in measurements of both the giant arc cross- 
section as well as the Einstein radius. In the previous sec- 
tions we have discussed the lensing properties of relaxed 
clusters so as not to be susceptible to the effects of massive 
substructures. We now consider the subsample of unrelaxed 
clusters to verify how their properties differ and how they 
would affect lensing predictions. We remind the reader that 
at zl = 0.25 there are 25 unrelaxed clusters and at zl = 0.5 
there are 20 unrelaxed clusters in our sample. Figure [14] 
shows the correlation between 07. 5 and 9e for clusters in the 
DM simulation comparing the relaxed and unrelaxed subsam- 
ple. The unrelaxed clusters in our sample would introduce a 
large scatter. We verifies that, as the threshold for the unre- 
laxedness parameter s max , is increased, the scatter increases 
in the relaxed sub-sample. This is primarily due to giant 
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in the CSF simulations are higher by ~ 20 per cent than 
for their collisionless counterparts. This boost increases to 
~ 50-70 per cent when we measure the Einstein radii at 
lower source redshifts. This is significantly larger than the 
boost for relaxed clusters. 

The introduction of AGN feedback once again reduces 
the strong lensing efficiency of unrelaxed clusters; the Ein- 
stein radii at z s = 2 for unrelaxed clusters in the AGN sim- 
ulations are only 5 per cent larger than those for their col- 
lisionless counterparts in the DM simulations. The difference 
in strong lensing efficiency is essentially the same for both 
relaxed and unrelaxed clusters. However, a mixed sample 
— with both relaxed and unrelaxed clusters — would have 
smaller Einstein radii on average, than a relaxed-only sam- 
ple. 



Figure 14. The relationship between the giant-arc cross section, 
0-7 5, and the Einstein radius, #e> f° r each of the 50 lines of sight 
analysed for each relaxed cluster in the DM simulation at zl = 0.5, 
for a source redshift of z B = 2. Results for relaxed clusters are 
shown in purple (solid line of best fit), while the unrelaxed clusters 
are shown in green (dotted line of best fit). The unrelaxed sub- 
sample exhibits a larger intrinsic scatter. 



arcs associated with substructure which induce a positive 
bias on the calculation of (77.5. The Pearson correlation co- 
efficient, r, for unrelaxed clusters in all simulations, ranges 
from 0.87 to 0.92 at z L = 0.25 and ranges from 0.82 to 0.90 
at zl = 0.5. 

The scatter above the <T7.5-#e lines of best fit is primar- 
ily because of substructures that are large enough to produce 
distinct caustics (see top row of Figure [TO")) ■ which induce a 
positive bias on the cross-section; as long as the caustics 
associated with projected substructures are well separated 
from the primary caustic, the Einstein radius measurement 
is not affected. On the other hand, if substructure are pro- 
jected near the clusters centre, their caustics merge and pro- 
duce highly elongated, critical curves (see bottom row of 
Figure [TO]), which artificially increases both measurements 
of strong lensing efficiency. However, the Einstein radius is 
affected more than the cross-section, thus down-scattering 
results for unrelaxed clusters with respect to the <T7.5-#e lines 
of best fit obtained for relaxed clusters. 

Figure [15] shows the cumulative probability distribution 
of Einstein radii for the unrelaxed sub-sample. Clusters iden- 
tified as unrelaxed are more likely to be merging systems, 
in which the main lens has a significantly lower mass than 
the total mass within the system. Hence, we find that typ- 
ical Einstein radii are only 50 - 70 per cent of the size of 
Einstein radii for the relaxed counterparts. 

We may now consider whether the effect of baryons on 
strong lensing properties is similar for the relaxed and unre- 
laxed sub-samples. Gas cooling and star formation produce 
clusters with a higher mass concentration; unrelaxed clusters 
have more substructures, which are also more compact than 
their collisionless counterparts (see, for example, the sub- 
structures visible in the cluster shown in Figure [5]). There- 
fore the unrelaxed clusters in the CSF simulations have larger 
primary caustics than their DM counterparts, but, compared 
to relaxed clusters, are also more likely to host secondary 
caustics that merge with the primary and boost the Ein- 
stein radius. Thus, the Einstein radii at z s — 2 for clusters 



4.4 Cluster Mass Profiles 

In the previous sections, we have seen that gas cooling, 
star formation and stellar and AGN feedback have varying 
and often counteracting effects on strong lensing efficiencies. 
Here, we look at the mass profiles of the cluster-haloes to 
trace the origin of the change in strong lensing properties. 
However, a detailed analysis of the mass profile fitting and 
the effect of baryons on the concentration-mass relation is 
beyond the scope of this paper. 

Previous analyses of non-radiative simulations in the 
literature have found that the gas pr ofiles exhibit cores (e.g. 
iRasia et all |2004| ; iRudd et al.1 120081 ). Figure [16] shows the 
spherically averaged differential density profile for the re- 
laxed clusters in two of our simulations: CSF and AGN. The 
profiles are shown for different components: dark matter, 
dark matter and gas, dark matter and gas and stars. The 
shaded regions mark the cluster-to-cluster scatter, while the 
error-bars mark the error of the mean; the small size of the 
error bars reflect the consistency in the profiles across the 
clusters. 

In keeping with previous analyses, we confirm that 
within r < 0.1R v ir gas cooling and star formation steepen 
the total mass profile with respect to DM-only simulations. 
In both CSF and AGN simulations, the baryonic component 
is dominated by stars within r < 0.1R v i r ; the gas compo- 
nent has a core and is distributed primarily beyond 10 per 
cent of the virial radius. The steepening of the inner profile 
is associated with the build-up of stellar mass; AGN heat- 
ing regulates star formation and therefore reduces the inner 
slope. We have mentioned earlier that comparisons between 
these two simulations must be interepreted carefully given 
that in the AGN simulations SN wind-speeds are increased 
while simultaneou s ly int roducing AGN feedback. Interest- 
ingly, iDuffv et al.l ((2J)f0j) have shown that in simulations 
without AGN feedback, the inner slope of the DM profile 
of clusters at z — actually steepens slightly when the SN 
wind-speed is reduced, but the stellar fraction within r < 
0.05R v ir is increased significantly, which we might expect to 
lead to an overall increase in the inner slope of the total mass 
profile. They find that the introduction of AGN feedback re- 
duces both the stellar fraction within r < 0.05R v ir and the 
inne r-slope of the DM p rofile. Our findings are in agreement 
with 
1§£ 



Duffy et al.l ( 2(ri0l) as well as several pre v ious studies 
Gnedin et al.ll2004l; IPuchwein et all 120051 : iMead et al.1 



l20I0l : iMartizzi et alJ^Oia T in wh ch the simulated clusters 
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Figure 15. The probability distribution for Einstein radii combining 50 lines of sight through each of the unrelaxed clusters in the 
sample. The source redshift is z s = 2. On the left panel, we show the results for clusters at z^ = 0.25; on the right panel, we show the 
results for clusters at redshift z^ = 0.5 




Figure 16. The differential density profile of dark matter (purple); dark matter and gas (green); dark matter, gas and stars (yellow). 
Clusters included in analysis are the relaxed sub-sample at zl = 0.5 in the CSF simulation (left panel) and the AGN simulation (right 
panel). The solid lines represent the average over all clusters, the shaded regions (in the respective colours) are between the 16th and 
84th percentiles among the clusters, while the error-bars correspond to the estimated error of the mean. The profile for each clusters has 
been normalised to the average mass density within the virial radius, and scaled with the square of the radius in order to reduce the 
dynamic range 



are analysed at z < 0.3. We included in our analysis also 
clusters identified at higher redshift, z = 0.5, and found no 
strong redshift dependence of the results. 

In order to see how the re-distribution of baryons affects 
the strong lensing, it is helpful to analyse the cumulative 
mass within different radii. In Figure[l7]we show the cumula- 
tive mass profiles for the clusters at zl — 0.25 and zl = 0.5, 
respectively. The smallest radius shown corresponds to 1 per 
cent of the virial radius, which corresponds to 16 h~ kpc 
for the smallest virial radius of all the clusters, or at least 3 
times the softening length. 

Our non-radiative NR simulations generate clusters with 
slightly more mass within r < 0.1R v ir than their DM coun- 
terparts; however, this is too mild to produce significantly 



strong er le nses, in agreement with the findings o f R ozo et al.l 
|200cf ) and lPuchwein et all (120051 ). The CSF simulations pro- 
duce clusters that con tain significantly more mass within 
r < O.lRvir (see also iGnedin et alj |2004| ; iPuchwein et al.l 
120051 . for comparable results on steepening of density profiles 
in radiative simulations). However, these simulations suffer 
the 'overcooling' problem — on average, 22% of the baryonic 
mass within rsoo in CSF clusters is in the form of stars — 
thus implying that the corresponding boost in strong lensing 
efficiency ought to be unrealistic. A more detailed compari- 
son between simulation and observation results on the stellar 
fraction will be presented by Planelles et al. (2012, in prep). 

AGN feedback reduces the total mass within 0.1 r v h- for 
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Figure 17. The cumulative mass profile for the relaxed clusters at zjj = 0.25 {left panel) and z^ = 0.5 (right panel). The lines reflect 
the average profile stacking all the clusters while the shaded regions are between the 16th and 84th percentiles among the clusters. The 
profile for each cluster has been normalised to the virial mass, and scaled with the radial distance in order to reduce the dynamic range 
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Figure 18. The Einstein radius relative to the cluster radius as a function of cluster mass. For each cluster, the median of the 50 analysed 
lines of sight are plotted with error bars reflecting the 16th and 84th percentiles over the distribution of lines of sight. On the left panel, 
we show relaxed clusters at zt, = 0.25, while relaxed clusters at zl = 0.5 are on the right panel. 



clusters at both redshifts 0. If we compare the DM and AGN 
simulations, we see that including AGN feedback results in 
clusters that contain more mass within r < 0.03R v ir relative 
to their collision-less counterparts. However this does not 
translate into larger strong lensing efficiencies. The reason 
for this is seen in Figure [TH] which shows the fraction of the 
virial radius encompassed by the Einstein radius at the clus- 
ter redshift, as a function of cluster mass. Typical Einstein 
radii are larger than the radius within which enclosed mass 
is affected by baryonic physics. So while clusters in the AGN 
simulations have a different structure with respect to their 



It is informative to consider AGN feedback prescriptions im- 
plemented in Adap tive Mesh Refinement (AMR) codes, such as 
those presented by iTevssier et al. | fcOllT ). They report a stellar 
mass fraction of /* fa 0.01 which increases with resolution to 
/* ~ 0.02. Although the concept of mass resolution is not di- 
rectly comparable between SPH and AMR codes, their results do 
caution that AGN feedback may be sensitive to resolution. 



DM counterparts, this difference does not produce stronger 
lenses, for sources at z s = 2. If we consider the results for 
sources at z s = 1, as shown in Table 1, the typical Einstein 
radii are smaller than for sources at z$ — 2. As a conse- 
quence, there is a higher significance to the probability that 
clusters in the CSF simulations are stronger lenses than their 
DM counterparts. Still, the strong lensing properties of clus- 
ters in the AGN simulations are still fairly consistent with the 
DM counterparts. 

We defer a more detailed comparison with observations, 
extending the simulated cluster sample to a wider range of 
masses, to a future study. For the purpose of the present 
study, it is sufficient to say that the suppression of star for- 
mation due to gas heating by AGN feedback in the AGN sim- 
ulations is responsible for the decrease in total mass in the 
cluster core, and subsequently, the decrease in the strong 
lensing efficiency of these clusters compared to their coun- 
terparts in the CSF simulations. 
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5 CONCLUSIONS 

In the present study, we have analysed the strong lensing ef- 
ficiency of about 15 relaxed simulated clusters at zl = 0.25 
and zl = 0.5. The clusters were re-simulated under a num- 
ber of different physical models: dark matter only (DM); 
non-radiative gas (NR.); cooling, star- formation, SN feedback 
(CSF); and additional AGN feedback (AGN). Our main find- 
ings are as follows. 

• There are no significant differences in the strong lens- 
ing properties and density profiles of clusters modelled with 
dark matter only and those that include a non-radiative gas 
component. 

• Gas cooling and star formation have the effect of in- 
creasing the giant arc counts and Einstein radii even in the 
presence of stellar feedback, which can push gas out of galax- 
ies and reduce star formation. 

• When AGN feedback is included in the simulations, the 
strong lensing efficiency of the clusters are significantly re- 
duced. Einstein radii at z B = 2 are only 5 per cent larger than 
those for the DM counterparts, a boost that is increased to 
10-20 per cent for lower source redshifts. There is a smaller 
orientation dependence, reflecting a more spherical shape 
due to gas cooling, even in the presence of AGN feedback. 
Small Einstein radii (8 < 3") are less likely to be produced. 

• The median Einstein radius is strongly correlated with 
the giant arc cross section for the relaxed cluster sub-sample. 
Unrelaxed clusters are especially problematic for the calcu- 
lation of strong lensing efficiency using the production of 
giant lensing arcs as a proxy; they also have smaller strong 
lensing efficiencies. 

• There exists a scaling relation between strong lensing 
efficiencies and mass, particularly for higher overdensities 
(A = 2500) . As we move to lower overdensities, the normal- 
isation for the each simulation differs primarily due to the 
redistribution of baryons near the cluster core. 

• We search for an explanation for these findings by 
analysing the mass profiles of our cluster sample; the dis- 
tribution of mass in the centres of clusters are significantly 
influenced by baryonic physics. AGN feedback reduces the 
stellar mass in the core of otherwise 'overcooled' clusters. 
Nevertheless, within r/R V ir 0.03 the total mass is still 
significantly larger than that of the DM counterparts. How- 
ever, this does not translate to higher lensing efficiencies for 
z 3 — 2, because 0e > 0.03R v ir for most clusters in our chosen 
mass range. 

Where there is overlap in cluster masses a nd redshifts, 
our re sults are similar t o tho se obtained by iDuffv et al.l 
120101 ) and iMead et all (|201Ch with regards to the effect 
of AGN feedback on density profiles and giant-arc cross- 
sections. The AGN feedback pres cription used in each work 
are all inspired by the model of ISpringel et all (|2005l ) but 
with somewhat different implementations and choice of the 
relevant model parameters. Therefore, it is reassuring that 
the strong lens ing predictions are n ot sensitive to these de- 
tails. However iTevssier et al.l l|201ll ). who include mechan- 
ical AGN feedback associated with jets into an Eulerian 
AMR code, have noted that stellar fractions are sensitive 
to implementation. An open question is then whether dif- 
ferent implementations of AGN feedback in different codes 
produce comparable effects on the inner density profiles of 



massive halos and, therefore, on the strong lensing proper- 
ties of galaxy clusters. 

As for the comparison with observational data, avail- 
able results on strong lensing arc statistics for samples of 
nearby and dist ant clusters h a ve be en recently carried out, 
for instance, bv lHoresh et al.l (|2010T l. Furthermore, the on- 
going Cl uster Lensing And Su pernova survey with Hubble 
(CLASH; IPostman et aHl2012T l project, thanks to the supe- 
rior sensitivity of deep HST imaging, is providing unprece- 
dented detail in the internal mass distribution of galaxy clus- 
ters through strong lensing studies. These results represent 
an excellent test ground for the comparison with simulation 
results, like those presented in this paper. Self-consistent 
comparisons with data are deferred for future study, as this 
will require more careful selection of clusters, with a com- 
parable mix of relaxed and unrelaxed objects in the obser- 
vational and simulated samples, as well as co nsistent meth- 
ods of measuring strong lensing properties. As lWvithe et al.l 
l|200lf ) have shown, a reduction in the lensing efficiency due 
to a decreased mass concentration may be partially compen- 
sated by a higher magnification for individual images. There- 
fore, a comparison of giant arc counts will require not only a 
prediction for the strong lensing cross-section, but also the 
magnification of images that are potential giant arcs, in or- 
der to determine if they would be detected in flux-limited 
imaging; this is particularly important for images that lie 
near the flux limit. 

Approaches to resolving the arc statistics problem have 
generally involved questioning three aspects of the cosmo- 
logical test: the assumed source redshift distribution, the 
realism of simulations, and the self-consistency in selection 
criteria for simulated and observed clusters. We prefer the 
Einstein radius as a proxy for strong lensing, since it allows 
one to avoid uncertainties with describing realistic source 
populations. None of the simulated models are tuned to re- 
produce all observable properties of clusters, but the AGN set 
of simulations is the most realistic one. Our study shows that 
AGN feedback is a game-changer; it's hard to explain a pos- 
sible discrepancy between observed and simulated lensing 
efficiency by resorting to the effect of baryons. Therefore, it 
will be crucial to make self-consistent comparisons between 
simulations and observations at a range of redshifts to un- 
derstand if lensing efficiency truly presents a challenge for 
ACDM. 
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